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^ ! Abstract 

We present a study of Bianchi class A tilted cosmological models admitting a proper nomothetic vector 
field together with the restrictions, both at the geometrical and dynamical level, imposed by the existence 
\Q • of the simply transitive similarity group. The general solution of the symmetry equations and the form 

of the nomothetic vector field are given in terms of a set of arbitrary integration constants. We apply 
the geometrical results for tilted perfect fluids sources and give the general Bianchi II self-similar solution 
and the form of the similarity vector field. In addition we show that self-similar perfect fluid Bianchi VIIo 
models and irrotational Bianchi VIo models do not exist. 
PACS: 4.20.-q, 4.20.Ha 
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The simplest generalisation of the Priedmann-Lemaitre universes are the Spatially Homogeneous (SH) models 
which admit a local group G r (r = 3, 4) of isometries acting on 3-dimensional spacelike hyper surf aces. From 



the physical point of view they serve as the basis in the study of anisotropics at early times of the universe 
bTj. evolution and play a key role in the understanding of the underlying geometrical/physical properties and 
j> | insights of more general cosmological solutions of the Field Equations (FE) : 

•i-H . 

Rab = T ab - —g a b- (1-1) 

_ _ I 

An extensively studied subclass of SH models are the Bianchi models where r = 3 which can be divided 
further into two major subclasses: orthogonal Bianchi cosmologies where the fluid four- velocity n a is orthog- 
onal to the group orbits [1] leading to the non-tilted models and Bianchi cosmologies where the normalised 
fluid velocity u a is not aligned with n a representing tilted models [2] . 

Although in SH models, the FE are reduced to a system of ordinary differential equations, not many exact 
solutions are known, especially in the case r = 3. This has initiated the examination of their dynamics using 
more sophisticated and predictive methods: the FE are reformulated as an autonomous system and using 
the methods of the theory of dynamical systems, the evolution of a specific model is studied in the so called 
state space which represents the set of all the physical states (at some instant of time) of the corresponding 
model [3, 4]. The main idea is to identify the evolution of a specific model with an orbit (solution curve) in 
state space and examine its geometrical-with direct physical interpretation-properties. One of the striking 
results of this approach is that transitively self-similar SH models are important in describing the asymptotic 
(i.e. at early and late times) states of the evolution for more general SH models. In particular, self-similar 
SH models arise as equilibrium points which determine various stable and unstable invariant submanifolds 
of the state space of the evolution equations. The determination of these invariant submanifolds allows to 
gain deeper insight into the intermediate evolution of the SH cosmological models. Furthermore there was 
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a strong evidence that any SH model could be approximated by self-similar SH models in the asymptotic 
regimes i.e. near the initial cosmological singularity and at late times. However it has been proved recently 
[5, 6] that Bianchi VIII and IX models are not (asymptotically) self-similar providing a counterexample in 
this conjecture. Nevertheless Bianchi type IX cosmologies are successively approximated (at the asymptotic 
regimes) by an infinite sequence of self-similar models (Kasner vacuum solutions) reinforce the belief that 
self-similar models have a significant role in the structure of the SH dynamical state space. 

It becomes clear that self-similar Bianchi models are of primary cosmological importance. We will focus 
attention on Bianchi models (r = 3) admitting a proper homothetic vector field (HVF) H acting simply 
transitively on space-time: 

£H9ab = tygab (1-2) 

where tp =const. is the homothetic factor which essentially represents the (constant) scale transformation 
of the geometrical and dynamical variables. 

The main goal of this paper is to provide the complete set of self-similar Bianchi models of class A 
(in particular for the "missing" types II,VIo,VIIo) irrespective the matter source of the gravitational field. 
In the case of non-vacuum fluid models the four-velocity subjects only to the assumption that is mapped 
conformally by the 4-dimensional similarity group. We utilize the purely geometrical results and study the 
case of tilted perfect fluid models. Due to the intrinsic complexity of the tilted models, up to date, not all the 
self-similar models have been determined adding to the difficulty of qualitatively analysing Bianchi models. 
Therefore the self-similar metrics presenting in this paper can be used to determine the general self-similar 
tilted perfect fluid Bianchi models of class A (note that in the vacuum and non-tilted case considerable 
information are available [7]). 

An outline of the paper is as follows: section 2 is devoted to a brief presentation of some of the basic 
properties of Bianchi models using the metric approach in which the basic variables are the frame components 
of the metric w.r.t. the canonical 1-forms [8, 9]. We then introduce a 1+3 decomposition of the tilted fluid 
velocity u a and show that the frame components of u a satisfy restrictions imposed by the structure of the 
isometry algebra of the SH models. In section 3 we confine our study to Bianchi class A models and present 
a set of constraints coming from the existence of a proper HVF. These constraints are used to determine 
explicitly the frame components of the metric, the fluid velocity and the corresponding similarity vector 
fields. In section 4 we discuss the physical implications of these, purely geometrical, results and examine 
the existence of tilted perfect fluid models. For Bianchi II models we regain the general self-similar tilted 
perfect fluid solution and prove that there are no irrotational self-similar tilted Bianchi VIq models. We also 
show that there no self-similar tilted perfect fluid Bianchi VIIq models. We conclude in section 5 discussing 
our results and further avenues of research. 

2 Preliminaries 

In this section we briefly review some basic features of the Bianchi models regarding the effects, both at 
the level of Geometry and Physics, which arise from the existence of the 3-parameter group of isometries. 
Since the theory is well known (see e.g. [3] and [10] and references cited therein) we present hear only 
those results which will be useful in the subsequent sections. Following Ellis et al. [3] we use the so called 
metric approach in which the basic variables are the frame components of the metric with respect to the set 
of canonical 1-forms. Although alternative techniques provide powerful unified formalisms of qualitatively 
analysing Bianchi models, the metric approach has the advantage of dealing with coordinate components 
hence we can determine the explicit (local) form of the self-similar metrics and consequently to study the 
resulting models directly. Throughout the paper the index convention is such that latin indices a, b, c, ...are 
tensor components and take the values 0,1,2,3 whereas Greek indices a, (3, 7,... = 1,2,3 denotes frame 
components (non-tensorial) of the corresponding quantities. A semicolon (;) denotes covariant differentiation 
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and a comma (,) denotes partial differentiation w.r.t. following index coordinate. 



2.1 Geometry 

Bianchi models are characterized by the existence of a G3 Lie group of isometries acting on 3-dimensional 
spacelike orbits C with generators the Killing vector fields (KVFs) X a tangent to the group orbits. The 
spacelike hypersurfaces C and the associated group G3 of isometries uniquely determine a unit timelike 
congruence n a (n a n a = — 1) normal to the spatial foliations C which satisfies: 

1 9 

n{a-,b] = = n a . b n b 44> -C n g ab = n a . b = a ab + -h ab (2.1) 

where <r ab ,6 = n a - b g ab ,h ab = g ab + n a n b are the trace-free symmetric shear tensor, the expansion and 
the projection tensor respectively associated with the timelike congruence n a [11]. It turns out that the 
(irrotational and geodesic) timelike congruence n a is invariant under the action of the group G3 i.e. the 
KVFs commute with n a : 



[X Q ,n]=0. (2.2) 

Therefore local coordinates {t, x a } can be chosen such that n a — — t a and the spacelike foliations C coincide 
with the hypersurfaces of constant t [1]. 

Introducing a set of invariant basis of vector fields e a and its dual to 13 such that: 

[X Q , e/3 ]=0 [X a ,u/]=0 (2.3) 

the full four dimensional metric can be written: 

ds 2 = -dt 2 + g al3 (t)Lj a LjP (2.4) 

and we have ensured the constancy of the 3- metric g a /3 on each orbit C. 

In addition the invariant basis e a and its dual u> a satisfy the corresponding commutation relations: 

[e a , e/3 ]=-C^e 7 , <kP = \cl^«N^ (2.5) 

where d stands for the usual exterior derivative of 1-forms, A denotes the usual exterior product 1 and 
C^p = —Cj a are the structure constants of the group. 

2.2 Physics 

As we have mentioned in the Introduction, attention is focused on perfect fluid sources which can be described 
in terms of the normalised four- velocity of the fluid u a (u a u a = —1): 

Tab = ifi + P)u a u b + pg a b (2.6) 

where fl,p are the energy density and the pressure measured by the observers comoving with the fluid 
velocity u a . Using the FE we may express the Ricci tensor in a similar form: 

Kb = (A + P)u a u b + — - — gab- (2.7) 



1 Note that the sign in (2.5) depends on the choice of the initial condition of the differential equations (2.3). 
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A straightforward consequence of the spatial homogeneity is that there is only one essential dynamical 
coordinate t. This follows from the invariance of the matter tensor under the action of the G3 [16]: 



Cx a T ab = £x a R ab = 0, [X Q ,u] = 0,X a (/i) = X«(p) = 0. (2.8) 

We note that the last two equations imply that perfect fluid models necessarily have a barotropic equation 
of state p = p(fji). Furthermore and for the purposes of the present paper, it is convenient to decompose u a 
parallel and perpendicular to n a as follows: 



u = Tn + YA a u a = Tn+TB a e a 



(2.9) 



where B a (t) , A a (t) are the frame components of the four- velocity u a and T is a smooth function of the time 
coordinate t satisfying the constraint: 



r = (1 - B a A Q ) 



(2.10) 



and B a A a < 1. 

The fluid velocity u a , although invariant under the action of the group of isometries, is not necessarily 
parallel to the normal timelike congruence n a , leading to tilted perfect fluid models i.e. the fluid flow lines 
are not orthogonal to the surfaces of spatial homogeneity (the quantity B a A a of the projection of u a is 
directly connected with the so called "tilt direction" [2]). In order to establish the relation between the 
dynamical quantities defined by the timelike vector fields n a ,u a , we consider the 1+3 decomposition of the 
matter tensor induced by the normal timelike vector field n a [11]: 



T ab = Hn a n h + ph ab + 2q (ya n b) + -n ah 



where 



(2.11) 



1, 



H = T ab n a n b , p = -T ab (n a n b + g ab ), q a = -h c a T cd n d , ir ab = h c a h d b T cd 

are the dynamical quantities relative to n a . 

Using equations (2.6) and (2.11) we obtain the relations [12]: 



1 



(h cd T cd )h ab 



(2.12) 



Ijl = jl + T 2 v\jl + p) , p = p+ -T 2 v 2 (il + p) 



q a = T 2 {fi+p)Y a , TT ab = T 2 (il+p) 



Y a Y b - hf 2 {g ab + n a n b ) 



where 



B a e a = A a co c 



y2 = yay = B a A 



(2.13) 

(2.14) 
(2.15) 



Equations (2.14) indicate that the energy flux vector field q a and the trace-free anisotropic pressure tensor 
ir ab are completely determined by the energy density and the pressure of the fluid and the spatial projection 
Y a of u a to the surfaces of spatial homogeneity. 

The system of equations is supplemented with the conservation equation: 



Tf = 

which for the case of the tilted perfect fluid source is reduced to: 



(2.16) 



V;aU a = -(fi+p)9 
haP;k = - {fl+p)Ua 



(2.17) 

(2.18) 
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and we have used the 1+3 decomposition of the first derivatives of the four- velocity [11]: 

9 ~ 

U a -b = &ab + ^h ab + Q ab - U [a ll b] (2.19) 

where a ab , 9 = u a - b g ab , Cl ab , u a are the corresponding kinematical quantities of the fluid velocity u a , namely, 
the trace-free shear tensor, expansion, rotation tensor and acceleration respectively and h ab is the projection 
tensor associated with u a . From (2.19) we obtain easily: 

U[a;b] = &ab ~ U [a U b] . (2.20) 

1 9 ~ 

^u9ab = U(a;b) = °ab + ^Kb ~ U( a U b y (2.21) 

Due to spatial homogeneity (third and fourth of equations (2.8)) all the (kinematical and dynamical) quan- 
tities are invariant under the G3 (which is equivalent with the vanishing of their Lie derivative with respect 
to the KVFs). 

The decomposition (2.9) implies: 

u a ;b = T yb n a +Tn a , b + r b A Q < + T (A Q ) >6 w a ° + rA Q < ;b . (2.22) 

We note that a similar relation holds if we apply the decomposition with respect to the invariant basis e a . 
The skew symmetric part of equation (2.22) gives: 

u [a . b] = - (TA a y n [b <] + TA a ^ a . b] (2.23) 

where a dot denotes differentiation w.r.t. t. 

Contracting equation (2.23) with e^e^ and using (2.5) we obtain: 

= ^C^A 7 . (2.24) 
In the case of a linear equation of state p = (7 — l)/i the conservation law (2.18) and equation (2.20) give: 

u [a . b] e a a e b p = (2.25) 

u [a u b] e a a e h p = 0. (2.26) 

Equation (2.24) relates the frame (spatial) components of ti[ a; 6] and Q, ab with the frame components A a of 
the four- velocity u a and the structure constants of the Lie group G3 (we note that a rather similar equation 
has been proved in [2] based on the Ellis's orthonormal tetrad formalism). Therefore a simple inspection 
of the group structure of the spatial homogeneous models, jointly with the FE, show that the kinematical 
and dynamical variables of the Bianchi models depend directly on the group structure. For example it is 
straightforward to prove that there are no tilted perfect fluid Bianchi I models whereas Bianchi types VIII 
and IX must necessarily have non-zero vorticity [2] . Indeed in type I the structure constants all vanish thus 
Qab = 44> U[ a;b ] = — ii[ a u b j which, by means of equation (2.23) and the evolution equation of Q ab implies 
ii a n a = 0. Then, provided that jx ]a 7^ 0, the conservation law (2.18) leads to u a = and the fluid is non 
tilted. In types VIII and IX, assuming that £l ab = equation (2.24) implies that A 7 = and the fluid is 
also non-tilted (and may be interpreted as anisotropic fluid). In the next section we draw attention only on 
Bianchi models of class A which possess a simply transitively self-similarity group. The case of class B will 
be studied in a forthcoming paper [13]. 
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3 Self-similar Bianchi models of class A 



A transitively self-similar Bianchi model admits a proper HVF H which together with the isometry group, 
form a simply transitive group of homotheties H4 acting on four dimensional orbits and has generators 
{H,X Q } spanned the associated Lie algebra denoted as H4. Our main concern in this section will be to 
find the general solution of the symmetry equations (1.2) in Bianchi class A models. However it is useful 
to present first the restrictions on the kinematical and dynamical variables of the Bianchi models induced 
from the existence of the HVF. For this reason we shall need the following well known result [14]: 
Proposition 1 

For every unit timelike vector field u a and a HVF H a = Mnu a + w a where w a u a = the following 
identity holds: 



L n u a 



(M H ). k u k + u k H k 



u k + 2uj ab H b + M] 



H 



(\nM n ),h k a 



(3.1) 



where ii a , u> a b, h k are the acceleration, the rotation and the projection tensor associated with the unit timelike 
vector field u a . 



3.1 Geometry 

A direct calculation shows that the commutator of a HVF with a KVF is always a KVF which implies that 
the structure constants C^ b of the homothetic Lie algebra TL4 satisfy C^ b = 0. Moreover the Jacobi identities 
for C^ b are equivalent to the relations: 

C] [a C L] = 0- (3-2) 
Following the considerations of the previous section, we decompose the HVF parallel and normal to n a : 

H = Hn + A a X a (3.3) 

where H and A a are smooth functions of the space-time manifold. 
From the commutator relations [X Q , H] = C2 4 X 7 we get: 

X Q (#) = , n(tf)=V (3.4) 



n(A') = , X a (AT) + ^Cj, = C2 4 (3.5) 

where we have used (3.1) and the fact that, since n is geodesic and irrotational, one has [H, n] = — n(H)n + 
[A a 'K a , n] = — ipn where n(H) = tp. Equations (3.4) and (3.5) determine the functions A 1 up to constants of 
integration. However some of them can be ignored by means of appropriate transformations which maintain 
the form of the metric and the KVFs. Therefore the HVF is easily determined, provided that C2 4 and 
C^g are given and subsequently, solving the symmetry equations (1.2), we determine the frame components 
g a p{t) of the metric. 



3.2 Physics 

Using the well known relations CuR a b = £-H.T a b = [15], it can be shown [16, 17] that the fluid flow lines are 
mapped conformally by the HVF and the equation of state is necessarily linear [18]. Moreover we deduce 
that in a self-similar SH perfect fluid model the "tilt angle" is constant and the projection Y a of the fluid 
velocity is also conformally mapped by the HVF: 

H(r) = o^r ;a = o , £ n Y a = -^Y a , £ H V a = VV a (3.6) 
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which by means of equations (2.22) and (2.23) imply: 

u a ;b = Tn a , b + T (A a ) j6 u£ + VA a u% b 
u [a -b] = -rA Q n [& <] + TA a uf a . b] . 
We may further observe that the frame components B a , A a satisfy: 

B a A a = const. 



(3.7) 
(3.8) 

(3.9) 



Equations (2.8) and (3.6) can be used to determine, up to a set arbitrary constants of integration, the form 
of the fluid velocity. 

3.3 Solution of the symmetry equations 

Solving the system of equations (1.2), (3.4), (3.5) and (3.6) we can find the explicit expressions of the self- 
similar metric (the frame components of the metric i.e. the functions g a /3{t) m (2-4)), the homothetic vector 
field and the tilted four-velocity u a for Bianchi class A models (we ignore the cases of Bianchi types I, VIII 
and IX since no tilted perfect fluid solutions exist). These results hold for a generic type of matter provided 
that, for fluid models, the four-velocity satisfies £h«" = —fpu a . 
Type II 

In local coordinates {t, x, y, z] and following the notations of [8], the KVs {X Q } and the canonical 1-forms 
{uJ a } are: 



Xi 



,X 2 = <9 Z ,X 3 = d x + zd y . 



(3.10) 

u 1 = dy — xdz ,uj 2 = dz , cj 3 = dx. (3-11) 

The non-vanishing structure constants are = 1 an d the Jacobi identities (3.2) imply that the remaining 
non- vanishing structure constants C^ 4 are [19, 20]: 



C\ A — a , C ; 



24 



r 2 



II- j ^24 

, <_/ 2 4 — 1 



C 14 — a , Cf 4 — 1 , Cf 4 — — 1 



1 ,Q 

•'34 " 



34 



,C 3 3 4 = a-l 
A 



34 



> C34 — a 



(a 2 < 4) 



(3.12) 
(3.13) 
(3.14) 



We consider subcases: 

Case Aj 
HVF 



H = iptdt + (a - 1 + b) d x + ayd y + zd z 



Fluid velocity 



Al = Vl t«>- a y+ , A 2 = t^-^ [v 2 + , A 3 = v^ +1 -^ 



(3.15) 



(3.16) 



Metric 



en 



t (2V>-a-l)/V 



Cl2 + 



fccn (In t) 



Cl3t (2V-2o+l)/^ 



if- 



i (2V>-a-l)/V 



C12 + 



t 2(V-l)/V 



C22 + 



6 2 



ben (In t) 



ciiQn ty 2bc r2 In t 
2 1 7, 



^(2V-a)/V 



C23 + 



6ci3(lnf) 



Ci3t (2V-2a+l)/V \ 

t (2ip-a)/ii> C23 + &siaj£j 

C 33t2(V-a+l)/V> 

(3.17) 



/ 
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where 6 = when a / 1. 

Case A2 
HVF~ 



H = 4>td t + xd x + X ^ Ay d y + (x + z) d, 



Fluid velocity 



Al = Vl t^~ 2 ^ , A 2 = v 2 t^- 1 ^ , A 3 = t^ 1 ^ (va - ^ In tj 



(3.18) 



(3.19) 



Metric 



9a/3 



Cnt 2(V-2)/V> Cl2 ^-3)/</> 



t 2{rl>-l)/rl> 
C33 



^23 " ^ 



(3.20) 



Case A3 
HVF 



H 



2-P2 



td t + (ax + z)d x +[ay + 



z 2 — x 2 \ 

— - — 1 d y - xd z 



(3.21) 



Fluid velocity 



Ai = v^- 1 , A 2 = t P2/2 
Metric 



Pi In t \ /pilnt 
V23 cos I — - — + V32 sin 



2-P2 



A 2 -t(A 2 ) )t (3.22) 



5a/3 



where: 



2-P2 



Cll t 2 fe-i) 
*(513) t - 2(p 2 - l)ffl 3 



521 



t (3p 2 -2)/2 J Ci;j CQS (pj^*) + gin ^ibit) 



2-P2 
2(2-p 2 ) 



i(523) jt -P2923 + 933 
t {933) t - 2 (p 2 - 1) 533 



531 \ 

532 

533 y 



(3.23) 



Pi 



,P2 



2^ - a 



■0 ^ 

533 = * P2 [C33 + C32 cos (pi In t) + c 23 sin (pi In t)] . 

We note that in all cases Ai,A 2 ,A3 the quantities v a ,c a f3,a,b are constants of integration. 
Type VIp 

The KVs {X a } and the dual basis {co a } are: 

Xi = d y , X 2 = d z , X 3 = d x + yd y - zd z . 
uj 1 = e~ x dy ,uj 2 = e x dz = dx. 



(3.24) 
(3.25) 



(3.26) 
(3.27) 
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In this case the non- vanishing structure constants of the isometry group C\ 3 = — Cf 3 = 1 and the Jacobi 
identities (3.2) imply that the remaining non-vanishing structure constants C^ 4 are: 



and the HVF assumes the form: 



C14 — a , Cf 4 — b 



H = iptdt + Dd x + aydy + bzd z . 



Furthermore the four-velocity and the metric are: 
Fluid velocity 



Ai = Vl t^ +D ' a ^ , A 2 = V2 t^- D ~ b ^ , A 3 = v 3 t 



Metric 



/ Cnt 2(^+D-o)/V Cl2t (2i>-a-b)/Tp Ci3t (2V+D-a)/^ \ 



9a/3 



(3.28) 
(3.29) 

(3.30) 
(3.31) 



Cv£t {2*l,-a-b)/i> C22t 2(V-D-b)/V C23t (2V>-D-b)/V 

V Cl3 t^+ D - a y^ C23 t^- D - b y^ " C33 t 2 j 

where v a , c a p, a, b, D are constants. We note that when the constants ci 3 , C23 are non- vanishing we may use 
the freedom y = ±y/ci 3 and z = ±z/c2 3 to set ci 3 = C23 = ±1. 
Type VIIq 



The KVs {X a } and the dual basis {uj a } are: 



Xi 



, X 2 = d z ,X 3 = d x - zd y + yd z . 



oj 1 = cos xdy + sin xdz ,uj 2 = — sin xdy + cos xdz , cj 3 = fix. 



(3.32) 
(3.33) 



constants C^ 4 are: 



The non- vanishing structure constants are Cf 3 = — C\ 3 = 1 and the remaining non- vanishing structure 

(3.34) 

(3.35) 



(7*4 — Cf 4 — a. 



The HVF has the form: 

H = iptdt + Dd x + aydy + azd z . 

The frame components of the four-velocity and the metric are: 
Fluid velocity 



Ai 



D\nt\ /Dint 
v\ cos I — ; — + V2 sin 



-v\ sin 



Dint 



A3 = v 3 t 



+ 772 COS 



Dint 



(3.36) 



Metric 



/ i 2(V.-a)/<0 
t 2(V>-a)/V 



C21 cos 
C21 sm 



/ 2D In A 



+ C12 sin 



/ 2D In A . i <_u 

v t J 2 



in(^)+ci 2 cos(^si) 



521 

i(2V-a)/V-[_ C31 sill (^*) + 
+Cl 3 COs(^' 



#31 
532 

C33* 2 



(3.37) 
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where v a , c a p, a, D are constants. 

Applying the translation x = x + A the form of the metric and the KVs remain invariant. Under this 
freedom in the choice of the x— coordinate we may set one of the constants 013,031, say C31, equal to zero 
i.e. C31 = 0. The other arbitrary constant, if not vanish, is fixed so that C13 = ±1 using the constant scale 
transformations y = ±y/ci3 and z = ±5/ci3. 



4 Tilted perfect fluid solutions 

Having found the general form of the Bianchi type II, VIo and VIIo metrics which admit a proper HVF, we 
may utilize these results to locate those space-times which have a tilted perfect fluid as their source. In this 
case the FE (1.1) reduced to a purely, although highly non-linear, algebraic system of equations. However 
using the restrictions on the frame components of the fluid velocity described in sections 2,3 we have been 
able to solve completely the system of equations for Bianchi type II and VIIo models and present some 
useful results concerning the Bianchi VIo models. 
The outline of the method is as follows: 

For each self-similar Bianchi metric we compute the Ricci tensor R a b and using the decomposition (2.7) 
we construct the tensor R a t, — jj}u a Ub— ^ 2 l^Qab which, together with (2.10), lead to 11 algebraic equations 
in 16 unknowns c Q/ g, a, 6, T, 7, jl, D,tp. We note that the number of the essential parameters needed for 
the general solution is reduced considerable for each Bianchi model as we shall see in the analysis presented 
below. 

Type II 

Tilted perfect fluid Bianchi type II models has necessarily zero rotation [2]. By virtue of equation (2.25) 
it turns out that Ai = 0. Furthermore, taking into consideration (2.9), the frame component A2 can be set 
so that A2 = 0. Therefore the general form of the tilted velocity in self-similar Bianchi II models is: 

u a = r (n a + ■ 



Analytical computations show that only in case Ai a perfect fluid solution exist which subsequently is 
the most general self-similar tilted perfect fluid Bianchi II solution. The metric is given by (3.17) with 
b = C13 = C23 = 0, C12 = 1 and the remaining non-vanishing constants are: 

_ 2p(3 - 4p) _ (1 ~ 4p) _ (4p-3)(5p-3) 

11 hp -I ' C22 2p(5p-3) ' C33 (4p-l)(5p-l) { ' 

where p is constant and we have set ip = ^3^, a = 3^7. 

The kinematical and dynamical quantities of this class of self-similar tilted perfect fluid models are: 

. = 2p(2p + 1) = J_ r _ y / 2(3-4p)(5p-3) 

M t 2 ' 7 2p+l' 2v/(5p-3)(3-4p)(9p-4) 1 ' ' 

where without loss of generality V3 = 1. 

Taking into account the overall signature of the space-time, the positivity of the quantities (4.3) show 
that these models are defined for p G (|, |) or p G (0, ^). The first set represents an exact Bianchi II perfect 
fluid solution where the parameter 7 lies in the interval [4/5, 10/11]. However this solution is not physically 
acceptable because the magnitude of the shear tensor a a ba ab and the measure of the acceleration ii a are 
both negative. In the second set of solutions the state parameter 7 G [10/7, 2] and has been found by Hewitt 
[21, 22]. 

Types VIp, VHp 

In contrast with the Bianchi type II models where significant information are available, little work has 
been done for tilted perfect fluid Bianchi VIo and VIIo models even in the special case of self-similar models, 
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due to additional degrees of freedom in the determination of the general solution. The problem is simplified 
with the help of equations (2.25) and (2.26) which, in conjunction with the FE and the generic form of the 
self-similar metrics (3.31) and (3.37), can be used to make some interesting comments concerning the known 
Bianchi VIo solutions. As far as we know the only known tilted perfect fluid solution is due to Rosquist and 
Jantzen [23, 24]. These Bianchi VIo models have non-zero vorticity and the parameter 7 £ (1.0411, 1.7169) 
(7 / 10/9). In the case of type VIIo we show that there are no non-vacuum self-similar solutions. 

From the above discussion one may wonder if there are non-rotating self-similar Bianchi VIo and VIIo 
tilted perfect fluid models at all. Although the isometry group G3 of Bianchi VIo and VIIo models, shares the 
common property of all the Bianchi models (except types VIII and IX) to admit an Abelian G2 subgroup, 
self-similar models are necessarily rotating. In fact we prove the following: 

Proposition 2 

There are no self-similar and irrotational tilted perfect fluid Bianchi type VIq and VIIq models. 
Proof 

Projecting equations C^ a u a = and Cuu a = tpu a with H and X a and substracting we obtain: 

2u [a . b] H a X b a = -^u k X k a + u k C^X k (4.4) 

which relates the projections of U[ a . b ] in the directions of the HVF and the KVs, with the structure of the 
homothetic algebra. Moreover equations (2.17) and (2.18) imply: 

u [a u b] H a X b a = ^l^± Uk X k a . (4.5) 

By means of equation (2.25), the condition Q ab = implies that the frame components A2 and A3 vanish 
therefore the KVs Xi,X2 are normal to the tilted velocity i.e. Xfu a = X^u a = 0. Since uy a .^ = —u^u^ 
from (4.4) and (4.5) it follows that 7 = 2 (we recall that C34 = 0) i.e. the only transitively self-similar 
Bianchi VIo and VIIo tilted perfect fluid models contain stiff fluid as its source. Then the FE (2.7), imply 
that the Ricci tensor has two eigenvectors Xf,X% with zero eigenvalue: 

R ab X\ = RabX b 2 = 0. (4.6) 

Hence the KVs Xi, X2 must be hypersurface orthogonal [25] and the metrics are necessarily diagonal. Taking 
the metrics (3.31) and (3.37) to be diagonal i.e. C12 = C21 = c 13 = C31 = C23 = a self-similar solution 
occurs only in type VIo (the type VIIo solution corresponds to the flat Friedmann-Lemaitre model). This 
solution is found to be unphysical (p < 0) in agreement with the general result that the only non-vacuum 
self-similar Bianchi model with an Abelian subgroup of isometries acting orthogonally transitively is given 
by the metric (3.17) and (4.2) [21]. 

In the case of rotating models our analysis shown that when 7 = 2 the resulting type VIo model is 
unphysical {fi < 0) hence there are no, physically acceptable, self-similar stiff fluid Bianchi VIo solutions. 
We have been unable to solve the resulting system of algebraic equations in full generality. However, 
although the family of solutions found in [24] belongs to the special classes satisfying the constraint = 0, 
we expect that this family represents the most general non-vacuum Bianchi VIo model. In fact setting 
a = (pi + P2) r 4> + 4^ — 6 and D = (2 — p<i) — b in (3.31) (pi,P2 are constants) it follows that the number of 
real degrees of freedom in a self-similar models is one (eleven equations and twelve unknowns). This issue 
is currently under investigation. 

Concerning the type VIIo models, the method described earlier in this section does not apply due to 
the complexity of the expression of the Ricci tensor. An alternative method to overcome the difficulties of 
computing the Ricci tensor, is the use of the Ricci identities for the timelike vector field n a combined with 
equation (2.26) and has as follows: 
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The first set of equations we will need is the constraint equations for the timelike congruence n a follow 
from the Ricci identities R a kbc nk = ^ n a[;cb] (equation (72) of [11]): 

h a c h b d (a cd ) b + q a = 0. (4.7) 
Contracting (4.7) with the KVs X® and using the first of (2.14) we obtain: 

(a%x£) a + r 7 fiuaXZ = 0. (4.8) 

In addition expressing equation (2.26) in terms of the KVs 2 we get: 

u [a u h] X%X\ = 0. (4.9) 

The set of equations (4.8) and (4.9) is used to examine the consistency of the self-similar metrics (3.37) with 
the assumption of a tilted perfect fluid source. Because the expressions of the KVs and tilted velocity are 
given by (3.32) and (3.36) respectively, we compute the quantity o a d X^ and the acceleration u a of the fluid. 
Replacing back to (4.8) and (4.9) we are left with a system of 6 algebraic equations whose solution gives the 
required answer. It is found that for every choice of the parameters, the energy density is always negative 
except in the case where the space-time reduced to the flat Robertson- Walker space-time. Hence there are 
no self-similar and tilted perfect fluid Bianchi type VIIq model except its flat RW counterpart. 



5 Discussion 

In this paper we have been able to give the explicit forms of the non-vacuum self-similar Bianchi models 
of class A regardless the matter content of the gravitational field. Therefore the results found in section 3, 
although purely geometrical, can be used to determine the general tilted perfect fluid self-similar solution 
since the only limited kinematical assumption we have made is the conformally mapping of the fluid velocity 
i.e. Cnu a = —ipu a . It follows that fluid models (non-tilted or tilted) with arbitrary energy-momentum 
tensor which admit a proper HVF, must included to the general models presented here. 

An important aspect of the analysis given in sections 3 and 4 (which is valid for the entire class of 
Bianchi models) is the direct dependency of the vorticity on the structure of the similarity group. This 
observation enables us to make some straightforward comments regarding the behaviour of the non-vacuum 
Bianchi models class B. In particular from equations (2.25), (4.4) and (4.5) we deduce that for a//irrotational 
Bianchi models, except type III, the fluid velocity frame components Ai = A2 = 0. Thus we expect that 
self-similar class B models subject to similar restrictions as in class A models. In fact it can be shown that 
self-similar and irrotational Bianchi IV and V models are only those for which 7 = 2 [13]. It turns out that 
the corresponding KVs Xi,X2 are also hypersurface orthogonal and the metric can be written in diagonal 
form in which case neither solution is physically acceptable hence there are no Bianchi IV and V self-similar 
tilted perfect fluid solutions with zero vorticity. 

Concerning the rotational class A models one must to solve the high non-linear algebraic equations in 
order to answer definitely whether or not a self-similar perfect fluid model exists. We have been able to work 
efficiently on the case 7 = 2 in Bianchi type VIo and show that this class of models, found by Wainwright 
et al. [26], do not admit a proper HVF. In addition based on the fact that there is only one essential 
parameter together with a rigorous, but not conclusive, study of the self-similar-metrics (3.31), we regard 
the assumption that the family of Rosquist-Jantzen solutions is the most general in type VIo, as highly 
plausible. 

2 We can always write the invariant vector fields as combination of the KVs in the form e a = Z^X-p where are time- 
independent smooth functions. 
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In the case of Bianchi type VIIo our analysis shows that self-similar models do not admit a perfect fluid 
interpretation. As we have mentioned in the Introduction types VIII and IX also do not admit a perfect 
fluid interpretation and it is well known that the corresponding non-tilted models share the same feature 
(see e.g. Table 1 of [7]). These results indicate that, in connection with the self-similarity breaking at late 
times in Bianchi VIIo non-tilted models [5, 27], it may be possible to relate the existence of a self-similar 
Bianchi model with the asymptotic self-similarity property of more general models. Consequently, it is 
an interesting and important task to determine under what conditions a full self- similarity breaking (i.e. 
in the asymptotic regimes) occurs and prove the conjecture that a general perfect fluid Bianchi model is 
asymptotically self-similar at late times or near to the initial singularity if and only if a self-similar perfect 
fluid model exists. 
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